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Abstract 

We have established a new convergent scheme to treat analytically nucleon-nucleon inter- 
actions from a chiral effective field theory. The Kaplan-Savage- Wise (KSW) amplitudes 
are resummed to fulfill the unitarity or right hand cut to all orders below pion production 
threshold. This is achieved by matching order by order in the KSW power counting the 
general expression of a partial wave with resummed unitarity cut, with the inverses of the 
KSW amplitudes. As a result, a new convergent and systematic KSW expansion is derived 
for an on-shell interacting kernel 7Z in terms of which the partial waves are computed. The 
agreement with data for the S-waves is fairly good up to laboratory energies around 350 MeV 
and clearly improves and reestablishes the phenomenological success of the KSW amplitudes 
when treated within this scheme. 
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1 Introduction 



Effective field theories are the standard method to deal with strong interactions in the non- 
perturbative regime. As paradigm we have SU(2) Chiral Perturbation Theory (CHPT) for pion 
physics 012 El 5 where a convergent power counting is established in terms of derivatives, insertions 
of quark mass matrix and external sources. This has also been applied to pion-nucleon interactions 
with baryon number, B, equal to 1 [I]. Its extension to nucleon-nucleon physics El HI EH El EH] 
is not straightforward due to the appearance of two new scales: the large scattering lengths of the 
S- waves and the large nucleon mass, M. The latter can be easily handled for B=l although this 
is no longer the case for B> 1 Large efforts have been devoted during the last years to end 
with a convergent effective field theory (EFT) for nucleon-nucleon interactions including pions, 
for a review see e.g. [11J. On the one hand, we have the original Weinberg's proposal [5 , with an 
undoubted phenomenological success In this scheme, the nucleon-nucleon potential is calcu- 
lated in a chiral expansion up to some definite order and then the Lippmann-Schwinger equations 
in partial waves are solved to determine numerically the physical amplitudes. Nevertheless, this 
approach suffers from inconsistencies in the power counting due to the appearance of divergences 
that are enhanced by powers of the large nucleon mass. As a result, their coefficients are much 
larger than those expected from 'naive power counting'. Despite this, such divergences are finally 
reabsorbed by counterterms whose chiral orders are established by assuming natural size following 
dimensional arguments. On the other hand, we have the Kaplan-Savage- Wise (KSW) effective 
field theory [T3] with its consistent power counting that directly applies to the physical amplitudes, 
like in standard CHPT pUEHH], and one finally has analytical nucleon-nucleon amplitudes. Inter- 
estingly, in the KSW power counting all the ultraviolet divergences appearing in loops are canceled 
by contact operators appearing at either the same or lower order in the expansion. This is not 
the case in the Weinberg's approach where one has cut-off dependence which gives an estimate of 
the size of higher order corrections. However, the convergence and the phenomenological success 
of the KSW effective theory is just restricted to a very narrow region close to the nucleon-nucleon 
threshold despite having included explicitly the pion fields [TH I15j. This was clear from the anal- 
ysis at NNLO in the KSW effective field theory performed in ref. ^1] where the NNLO departs 
from data before the NLO and badly diverges for center-of-mass three-momentum above ~ 100 
MeV. The reason for this bad convergence properties is the large contributions from the twice 
iterated pion exchange in the triplet channels. That is, pions cannot be treated perturbatively in 
all the nucleon-nucleon channels and in particular in the Sf — Df one. 

Hence, despite all the efforts made during the past years, the issue of deriving a consistent 
effective field theory for nucleon-nucleon interactions is still open. We perform in this paper a 
step forward by extending the range of convergence of the KSW amplitudes. In ref. Wein- 
berg proposed to calculate the nucleon-nucleon potential in a chiral expansion and then solve a 
Lippmann-Schwinger equation in order to treat properly the large enhancements, due to factors 
2M/p, from the reducible two nucleon diagrams which do not enter in the calculation of the po- 
tential, see also ref.jH]. From a S-matrix point of view the aim of the Weinberg's proposal can be 
recast so that one should resum the unitarity or right-hand cut, which is the responsible for all 
the unitarity bubbles with their large 2M/p factors. Nevertheless, solving a Lippmann-Schwinger 
equation is not the only way to accomplish this ^H] and another one, more appropriate for quan- 
tum field theory can be established. In this way, one avoids the non trivial problems associated 
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with the renormalization of a Lippmann-Schwinger equation with highly singular potentials at the 
origin as those coming from present effective field theories. 

We propose here a general scheme to resum the unitarity cut. This scheme, originally motivated 
as an application of the N/D method [HI HE], was already employed in meson-meson [THIHH] and 
meson-baryon fZU\ !2Tj systems. This method should not be confused with Pade resummations 
like those performed in the meson-meson or meson-baryon sectors [22j-* 2 Nevertheless, there 
are specific facts in the nucleon-nucleon scattering, associated with the largeness of the scattering 
lengths in the S- waves, that require a special treatment not present in any meson-meson or meson- 
nucleon system, as we explain below. In fact, a very reassuring feature of nucleon-nucleon physics 
is that non-perturbative effects manifests at very low center of mass three-momentum, p, as 
indicated by the presence of bound states and poles in unphysical sheets just below threshold. 
Hence, one has at his disposal the three-momentum p as a small parameter and this is of course 
at the basis of the old Effective Range Expansions (ERE) 23j. It is also worth stressing that all 
the formalism presented here is analytic. 

The paper is organized as follows. In secJ21the novel formalism contained in this work is settled. 
This follows the KSW power counting and the reader unfamiliar with such power counting is 
referred to the original literature, refs. |13| I14j. It is shown how the unitarity cut can be easily 
resummed in terms of a dispersion relation of the inverse of a partial wave amplitude. As a result 
an interacting kernel TZ arises that is determined by matching with the pure KSW amplitudes 
following the KSW power counting. Once this is done, we discuss in secJ3] the phenomenology 
and compare with data. Since all the expressions are analytic, we fix, order by order in the 
expansion, many of the counterterms in terms of the ERE parameters a s and r in the S-waves. 
The agreement with data for the S-waves and the mixing angle Sf — Df is quite remarkable in a 
broad energy range, p < 400 MeV (7] a 6 < 350 MeV), involving all the data points of the Nijmegen 
partial wave analysis [21] • Higher partial waves are also analyzed at the same order. However, only 
one pion exchange and the reducible part of twice iterated one pion exchange enter in these KSW 
amplitudes up to order p. The results indicate that higher order corrections should be included 
so as to perform a more complete analysis for the P and D partial waves, as in the Weinberg's 
scheme jTS]. We end with some conclusions. 

2 Formalism 

It seems that pions cannot be treated perturbatively. This is indeed the main reason for the failure 
of convergence of the KSW scheme, which considers pions perturbatively in the physical ampli- 
tudes, for center of mass three- momentum p > 100 MeV. This was clearly established in ref. [14]. 
as one can see when comparing the calculated mixing parameter e\ at NLO [IS] and NNLO [T4] . 
While for very low momentum the NNLO calculations agree better with data, unfortunately for 
momenta higher than 100 MeV the NNLO calculations badly diverge from experiment, much more 
than the NLO results, and no improvement is obtained despite pions being explicitly included. 
This should be compared with the phenomenological success of the Weinberg's scheme in jH] and 
particularly in where the chiral expansion seems to be under control, at least at the level 

# 2 Later one we will show that these Pade resummations are particular cases of our formalism. 
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of the phenomenology. In the latter scheme there are still the issues of avoiding cut-off depen- 
dence and the associated inconsistencies in the Weinberg's power counting, already referred in the 
introduction and originally established in refs. |Hl HH1 • 

In this paper we make use of the KSW amplitudes which are calculated within an effective field 
theory and are properly renormalized. There is nonetheless a clear difference between the KSW 
effective field theory and CHPT. While the latter is convergent in the SU(2) sector in an energy 
window that could be expected from the scales that are involved in the problem,* 3 the same does 
not occur with the former |14j . 

We want to include pions non-perturbatively but keeping at the same time the advantages of the 
KSW scheme that make it a true effective field theory. In order to accomplish this, we resum the 
right hand cut or unitarity cut so as to take care of the large 2M/p factors from the two nucleon 
reducible diagrams [3]. This is also performed when solving a Lippmann-Schwinger equation, as 
originally proposed in ref.jH], but we will do it in a different fashion considering unitarity and 
analyticity in the form of a dispersion relation of the inverse of a partial wave. 

Let us denote a partial wave by T L 2s+i ps+i , where L and L refer to the orbital angular mo- 
mentum of the initial and final state respectively, and S and J indicate the total spin and total 
angular momentum of the systems, in order. Then unitarity requires above the nucleon-nucleon 
threshold and below the NNtt one, which is around p ~ 280 MeV: 

ImT Lf+\Lf +1 = ^^ T Lf+\e 2 /+ 1 ( T e 2 J s + 1 ± 2 J s+1 Y > (2- 1 ) 

i 

where M is the nucleon mass. It is easy to derive from the previous equation: 

Im(T ( 2S+ iy).i = _^, (2 . 2) 

where T^ 2S+1 ^ J is a matrix whose matrix elements are those triplet partial waves that mix each 
other, e.g. for Sf — Df one has Tfl = T S 3 S 3, = Tfi = ^Sf,Df an d = T D 3 D 3. In the previous 
equation we denote by (T < - 25+1 ) J )~ 1 the inverse of the T^ 2S+1 ^ J matrix. If the partial waves do 
not mix then J"( 25 + 1 ) J is just a number equal to T L 2s+i L 2s+i. The imaginary part in eq. (J2.2|) is 
responsible for the unitarity cut. Thinking of a dispersion relation for the inverse of the amplitude 
this cut can be easily taken into account from eq. (|2.2|) . as we previously did in the meson- meson 
[T5] and meson-baryon [20J |^j sectors, and gives rise to the integral: 

, , 1 f°° Mp(s') 1 . . 

9i(s) = ~ / Hl ~ r^T ds . (2-3) 

vr J 4M2 47T s' — s + i0 + 

where s is the usual Mandelstam variable. This integral is divergent and requires a subtraction, 

, > M f . Ma(s) , l-a(s)\ 

gi (s) = — U + ip+ ^log , 2.4 

An \ i\ 1 + er(s) / 



with cr(s) = y 1 — ^Ml. The logarithm is purely real in the physical region and just gives rise 
to relativistic corrections. These are obtained from eq. (|2.3j) but with Mp(s') / 'Att replaced by 
# 3 Except in the scalar isoscalar channel. 
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the relativistic phase space M 2 p(s') /2n^J~s in the integrand. Since we want to match with the 
KSW amplitudes, that absorb all the relativistic corrections in the vertices, the non-relativistic 
phase space factor Mp/4-7r is used as explicit source of imaginary part in eq. (|2.4jl . The relativistic 
corrections from the logarithm of eq. (|2.4jl are in any case essentially negligible in the energy range 
that we consider. Once the unitarity cut is taken into account by the functions <?;(s), the full 
partial wave matrix or number T ( - 2S+1 " )J can be written as: 

T (2S+l)J = _^-l + ^-l (25) 

In this expression all the other possible cuts of a nucleon-nucleon partial wave, either due to the 
exchange of other particles (in our case we have the exchanges of pions) or because its helicity 
structure, are included in the interacting kernel TZ that we still must fix. The unitarity require- 
ments, resummation of the infinite set of reducible diagrams with two intermediate nucleons, are 
accomplished by g(s), which as stated is a diagonal matrix in the case of mixed partial waves, or 
just one function for the unmixed ones. For example, in the coupled partial waves Sf and Df, 
9n( s ) = 9i( s ) an d corresponds to the Sf channel and (722(5) = #2(3) and refers to the Df channel. 
Of course, g 12 {s) = g 2 i(s) = 0. 

We now specify TZ in the key expression eq. (|2.5|) . For that purpose we make use of the results 
of the KSW effective field theory for two nucleon systems and of its power counting, that we apply 
to TZ and g(s) in eq. (|2.5|) . The matching procedure with the KSW amplitudes can be done for 
any given order in the calculation of the KSW amplitudes and for any subtraction constants V{ in 
<7i(s), as we show explicitly below. 

An important point is to establish the chiral order of gi(s). This is a trivial task for the phase 
space and the logarithmic term in g(s), eq. (|2.4jl . since they can be expanded as a series in powers 
of p starting at first order. The only point we have to consider separately is the chiral order of 
z/j. For the S- waves Sq and Sf, the corresponding elastic KSW amplitudes start at order p -1 and 
hence their inverses begin at order p. Thus one can take V{ as order p as well, since this is the 
first order that appear in the inverses of the leading KSW partial waves. Would this be the case, 
then TZ would also start at order p and the matching with KSW amplitudes can be performed 
straightforwardly. However, if we continue along these lines there is no clear improvement with 
respect to the problematic of the KSW calculations at NNLO. If, on the other hand, we think of 
higher partial waves, e.g. P and D, then we realize that the KSW amplitudes start at order p° 
so that taking g and TZ to start as order p° is quite natural. If we take this option as well for 
the S-waves, and systematically derive the TZ matrix elements by matching with the KSW inverse 
amplitudes, we will see below that the improvement with respect to the pure KSW scheme is 
fairly remarkable. Indeed, in order to match with the inverses of the KSW S-wave amplitudes, 
that start at order p, we must cancel exactly the order p° contributions in eq. ([2.5|) stemming from 
the ones of TZ and g. This fine tuning reminds of the one usually advocated to explain the large 
scattering lengths in the S-waves channels [3]. 

In order to clarify further the previous discussion, we can easily see that considering v\ as a 
constant of order p° is a result when g(s) is calculated with a finite cut-off. Let us perform this 
illustrative exercise of more than academic importance since the finite three-momentum cut-off is 
the regularization employed in the Weinberg's scheme. For a given nucleon-nucleon channel i, we 
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MeV 


MeV 


MeV 


MeV 


MeV 


MeV 


100 


64 


500 


310 


900 


510 


200 


130 


600 


360 


1000 


550 


300 


190 


700 


410 


1100 


600 


400 


250 


800 


460 


1300 


670 



Table 1: Values for Vi from eq. (|2.10j) for different values of the three-momentum cut-off Q. 



can write: 



gt{s) = —iAM 2 



d A q 



(2vr) 4 (q 2 — M 2 + z0+)((P - q) 2 - M 2 + i0+) 



(2.6) 



where P is the total four momentum of the two nucleon system, P 2 = s, and the superscript 
c in g c (s)i indicates that is calculated with a three-momentum cut-off. After performing the q° 
integration we have: 



M 2 



w- 



M 2 



2tt 2 J m w 2 - s/4 - i0+ 



dw 



(2.7) 



where A = \/ Q 2 + M 2 with g the three-momentum cut-off. Finally we obtain the explicit result: 



m 2 I , A + g 



log 



2M 2 +A-y/5 



log 



2A 



(7 



2A+ v/i 



f2.8) 



Performing a non-relativistic expansion in the previous equation we have: 



M 2 
i^ 2 



A + g .7Tj9 1 — a(s) p 2 , 

2 log — + «— + o-(s) log H- + 0( At 

M M 1 + <x(s) V M 2 ' 



(2.9) 



Comparing with gi(s), eg. ()2.4j) . we finally have: 
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2M . A + g 

log 



7T 



M 



(2.10) 



which is a quantity of order p° in the KSW power counting [T3*l ITlj . In the KSW EFT Q is 
expected to be around 300 MeV, [T3] . In table [T] we show the values of u \ in MeV from eq. ()2.10|) 
for different values of Q. We will obtain later, directly from fits to data, similar values of v^. 

Before applying the previous scheme to the KSW amplitudes, let us consider the pedagogical 
example of the expansion in powers of x of f(x) = cotx = cosx/ sinx. For that we write: 



t{x)- 1 +6{x) 



(2.11) 
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such that t(x) = t + t\x + t 2 x 2 + 0(x 3 ) and 6(x) = z + Z\x + z 2 x 2 + 0(x 3 ). We see here that 
although f(x) starts at order x _1 we have considered the functions r and 9 to start at order x°. 
In order to fix ti in terms of the Zi (which are assumed to be known) and of the known expansion 
of cotx, is simpler to expand the inverse of f(x), then we obtain up to 0(x 2 ): 
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= x + 0{x 3 ) = -— + %x + %x 2 - l -\x l -zq- z x x - z 2 x 2 + 0{x 3 ) . (2.12) 
cos x to £q £q /^q 

It follows then: 



l 

{zi + 1) 



Z 2 (Zi + l) 2 , 
*2 = 4+ 3 • ( 2 - 13 ) 

It is obvious how to proceed for higher orders. This simple example also illustrates that if we 
want to calculate t(x) up to order x 1 one needs to know f(x) up to order x l ~ 2 since f(x) already 
starts at order x~ x . 

Let us now analyze carefully the elastic S\ channel following the previous scheme. After that 
we will present more briefly the analogous procedure in the Sf — Df coupled channel sector and 
for the P, D, F 2 2 and G% waves. 

2.1 Sq elastic partial wave 

The KSW partial wave, A^ sw , was calculated at NLO (order p°) in ref. [T3] and then at NNLO 
(order p) in ref. [Hj. Let us denote these partial waves by A_i(p), A (p) and Aj(p), where the 
subscript indicates the KSW order. As in the simple example of the cot x, if we take as input 
the KSW amplitudes up to order p then we will be able to calculate 1Z up to order p 3 . This 
unambiguously fixes the order one has to calculate in the KSW EFT so that TZ is obtained up to 
the required precision. Following the same notation as for the KSW amplitudes, let us write: 

K = R + Ri + R 2 + Rs + 0(p A ) , (2.14) 

and g(s) = Mis/Air + iMp/Air - p 2 /2ir 2 + 0(p 4 ). Then we we must match: 



1 / 1 \ Mo V (Al-A X A_, 



A* sw \A_J \A 2 ' V A 3 



0(p 4 ) , (2.15) 



with 



i , \ = _(1_ + Ml] + (?±- iME] + (f- + RoR2 ~ R2 ^ . ~ 2R ° R i R i + ^o^s 



n y \Rq 4tt J \R 2 Q Air J \2ir 2 R 3 J \ R A 

+ 0(p A ), (2.16) 
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where we have shown between brackets the different orders, from order p° up to order p 3 . As a 
result of the matching we can fix _R , Ri, R% and R3 in terms of u, A-i(p), A (p) and A\(p). 
We follow ref.jHj for expressing the KSW amplitudes, so that any of the previous amplitudes are 
scale independent at each order in the expansion. 
Taking into account that 

4?r 1 . 

A-i = "TV > 2.17 

M 7 + ip ' y 1 

with 7 a quantity of order p, we can write the following expressions for the R4: 



Ro 
Ri 



R 



2 



An 

~Mv ' 

477T 

4(2i^p 2 + 1 2 Mn) + u(^- i ) 2 A (p) 



1 Att Att 

R 3 = -^^^^TvA^-v^^f 

+ A^TMP)]) ■ (2.18) 

Hence working at NLO, 0(p°), in the KSW amplitudes [T3] we will have 1Z up to 0(p 2 ), as 
explained above, 

TZ NLO = R + R 1 + R 2 , (2.19) 
and matching with the ones at NNLO [13], 0(p), we calculate 1Z up to 0{p 3 ): 

n NNLO = Rq + Ri + R2 + R3 (2.20) 



The resulting 1Z is substituted in eq. (|2.5jl and in this way we calculate the partial waves at the 
different orders considered so forth. It is clear that the process above can be done so as to match 
with a KSW amplitude calculated at any order. In this way the precision of the resulting amplitude 
is increased order by order. 

Let us consider now in more detail the meaning of the kernel 1Z. Take first the simpler case of 
the pionless effective field theory, where only local interactions and even numbers of derivatives 
appear in the Lagrangians. For this case, the resulting partial waves have only the right hand cut 
and are free of crossed cuts due to the exchanges of other particles, being the pions the lightest 
ones. This is important for the nucleon-nucleon dynamics since at very low energies the pions can 
be treated as heavy particles and integrated out j2H] • Thus, only local operators and the unitarity 
cut remain. We follow the line of reasoning of ref.[2Z]; where the so important Castillejo-Dalitz- 
Dyson (CDD) poles were introduced. This reference was also used in ref. [IB] in the context of 
chiral Lagrangians to show the general structure of a meson-meson partial wave with only the 
unitarity cut. Here we show the main points of the reasoning, for further details the reader is 
referred to ref. JH]- The idea is that when the partial wave T L 2s+i L zs+i has a zero then its inverse 
has a pole. This pole in the inverse of the partial wave can be on the real axis or in other place 
of the complex plane of the physical sheet (Imp > 0). Special care is needed when these poles lie 
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on the real axis above threshold since then eq. (j2.1|) can only be applied between any pair of such 
poles. These technical details are given in refs. [23 EI] • In this way, the inverse of the partial wave 
is a meromorphic function in the cut s plane from threshold to infinity. The net answer from a 
dispersion relation of the inverse of the amplitude, by applying the Cauchy theorem of integration 
to the circuit made up by the circle at infinity deformed to engulf the unitarity cut, is: 

7: 



t l ?»w = - Ettt+H ' (2 - 21) 



where the poles present in the sum are the so called CDD poles [27,. The sum is what we have 
denoted before in eq. ()2.5j) by 1/11. Let us note that in the pionless effective field theory supplied 
with the PDS scheme, the Schrodinger equation can be solved straightforwardly and results |lHj : 



with A being the same partial wave as T L 2s+i L 2s+i. Comparing this equation with eq. (|2.21j) . one 
has: 

1 _ 7n 1 M _ 1 

where we have omitted in the previous equation the relativistic corrections in g from the logarith- 
mic term in eq. (|2.4|) . The C' m are given in terms of the C rn so as to reabsorb the constant term in 
the sum. One can always choose the s n and 7 n so that the previous equality holds.* 4 . This is why 
we have written finally 1/1Z in eq. fl2.5j) . instead of simply 1Z, because at the tree level, omitting 
the g function in eq. fj2.21j) . one has T L 2s+i L 2s+i = 1Z. Let us note in addition that —1Z can be 
identified with the renormalized potential if we rewrite T = — 1/ (1/1Z + g) as: 



T = -1Z-1ZgT, (2.24) 

like an ordinary Lippmann-Schwinger equation with an on-shell potential V = —1Z. Let us remind, 
as stated above, that g is just the unitarity bubble iterated when solving a Lippmann-Schwinger 
equation. 

When the pions are included in the formalism, the interacting kernel 1Z contains apart from a 
finite sum of local terms, like those of eq. (j2.2Hj) . other contributions coming from the exchange of 
pions which give rise to crossed cuts. All these contributions are included perturbatively when 1Z 
is fixed by matching the expansion of the inverse of eq. (j2.5j) with the expansions of the inverses of 
the KSW partial waves. Let us stress that the entire formalism is algebraic since the interacting 
kernel 1Z is on-shell. This notorious simplicity of our approach compared to that of refs. El El 
is a result of resumming the unitarity bubbles by making use of analyticity and unitarity that 
only involve on-shell amplitudes. One needs to realize that for resumming the large factors 2M/p 
associated with two-nucleon intermediate states, whose necessity was stressed in ref. jS], there are 

# 4 For example, prove the above equality by considering first that the sum on the right hand side just contains 
two terms. Then, for the general case, it is easily seen, by recurrence to the case with one term less in the sum, 
that the equality holds. Note that when s n — > oo the corresponding CDD pole just gives rise to a constant 
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several possibilities. One way is to solve the Lippmann-Schwinger equation OIHIEIj but there is 
another standard method based on the separation of cuts resulting from S-matrix theory, this is 
the N/D method [Tj\ that we follow here. 



2.2 Sf—Df coupled partial waves 

The resulting expressions to match between are the expansions of the matrix elements of the 
inverse of the KSW matrix of partial waves for the Sf — D\ sector, A^^: 

^12,0 ~~ ^11,0^22,0 

^11,-1^22,0 

^12,0 + 2y4 llj _ 1 A 12)0 A 1 2 ; lA 2 2,0 + (^11,0 — ^H,-l^ll,l)^22,0 

43 42 

^11, -1^22,0 

^12,0(2^11,0^22,0 + All, -1^22,1 

43 72 
/1 11,-1 /1 2 



+ 



(A 



+ 0(p 4 

22,0 

^12 \ ( ^12,0 + ^11,-1^12,1^22,0 ~~ ^12,o(Al,0^22,0 + Al,-1^22,l' 



^12,0 — ^11,-1^22,1 \ . nf 2 \ 



2 

22,0 



and the ones that result from the expansion of eq. (|2.5|) . In this case 1Z is a 2 x 2 symmetric matrix: 

* = f ^ ) , (2-26) 

V -K12 -K22 / 

and g(s) = diagonal(gi(s) , g2{s)) with its associated subtraction constants V\ and z/ 2 . Like in the 
S\ channel we will take g t = O(p ) and i?n = #u j0 + #u,i + R llt2 + i?n j3 + C(p 4 ). In the KSW 
scheme at the leading order the Sf and .Df are uncoupled and the mixing starts at NLO, one order 
higher. Thus, we take i?i 2 = #12,1 + #12,2 + 0(p 3 ), starting one order higher than i? n . Finally, 
since the Df partial wave starts at order p°, and is free of unnatural scattering lengths, we take 
then R 2 2 = -^22,0 + -^22,1 + 0{p 2 ). The expansion of the inverse of eq. (|2.5|) is now straightforward 
and one obtains: 

/rrVm-l ( M.V\ 1 \ ( R\\\ Mp\ ( P 2 -Rll,0-Rll, 2-^22,0 — -Rll 1-^22,0 — -^11,0-^12 1 

( T hi = + ^ — + 772 1 



4vr ' Ru fi J ' V-Rii.0 4?r / \2vr 2 #11,0-^22,0 

-^11,1-^22,0 + 2i? llj0 -Rii,i-R22,o( — -R?2,i + -Rn, 2-^22,0) 



d4 d2 
- fl ll,0- rL 22,0 

#11 o( _ 2-Rl2, 1-^12,2-^22,0 + #11,3-^22 + #12 1#22, l) 



R4 d2 
- rL ll,0- rt 22,0 



+ C(p 4 



-,31\-1 / #12,1 \ / #11, 1#12,1#22,0 ~ #11,0#12,2#22,0 + #11,0#12,1#22, 



'■" .#11,0#22,0/ V #11,0#22,0 
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We can then easily solve for the Rij,k an d we obtain: 

47T 



Rn,o — 

RllA = ~M^' 

Rn, = M2 ,3 (47r I Mu2A22 o) {^ 2 (7 + i P ) 2 A* 2fi - 4(2u lP 2 + 7 2 Mir)(4ir + Mu 2 A 22fl ) 

- M 2 v 1 (- i + ip) 2 A llfl (An + Mv 2 A 22fi )} , 

i?11 ' 3 = 4Mu 2 7iA 2 22 {" M(7 + + ^) Ao " 8nA 12fi A 12A A 22t0 + (M( 7 + ip)A\ lfl 

+ 47rAn,i)^22,o " 2A 2 2>0 (M( 7 + ip)A llfi A 22fl - 2nA 22A )) + 1 ( 

^ (47T + M ^2 ^122,0 J 

x 87r(2 7 3 7r^ 2 (47r + Mz/ 2 A 22i0 ) 2 - -^( 7 A 22i0 (47r + M l / 2 A 22 , ) 2 (-MV 1 ( 7 + ip) 2 A 2 20 

+ {4{2v lP 2 + 7 2 Mvr) + MV( 7 + ip) 2 A lli0 )A 22i0 )) - + ipf A 12i0 (-Mv 1 (>y + ip)A\ 2fl 

x (2vr + Mz/ 2 A 22i0 ) + 4z/i7rAi 2i iA 22i o(47r + Mz/ 2 A 22i0 ) + Ai 2i0 (M(4 7 z/ 2 7r - 2iuipn 
+ Mi/ 1 i/ 2 ( 7 + ip)A llfl )A 2 22 Q - %vi-k 2 A 22)1 + 47rA 22i0 (4 7 7r + Mv^ + ip)^n, 

- Mu 1 u 2 A 2 2,i)))))} , 

4(7 + ip)7rAi 2;0 
12,1 v^Air + Mv 2 A 22fi ) ' 

#12,2 = 271 T a# 4 ^{(7 + ^)(-M 2 z/ 1 z/ 2 ( 7 + ^)^ +4z/ 1 7tA 12i1 (47T + Mz/ 2 A 22i0 ) 

^47T + 714 Z/ 2 /i 22i o J 

+ ^12,0 (M1/1 (7 + ip)Aii )0 (47r + Mi/ 2 A 22i0 ) + 47r(4 7 7T + M( 7 z/ 2 - iz/ip)A 22j0 - Mi/ii/ 2 A 22) i)))} , 
47tA 22;0 
22 '° 4tt + Mz/ 2 A 22i0 ' 

4tt(M( 7 + tp)A 2 2fi - «MML,o + 4ttA 22i1 ) 
i?22 ' 1 " (4 7 r + M, 2 A 22 , ) 2 • (2 ' 28) 

Similarly as in the Sq case, working at NLO implies: 

Rn L ° — Rn,o + Rn,i + Rn,2 , 

piVLO _ p 
-"■12 ~~ -n-12,1 ) 

-^22 L ° = -^22,0 , (2.29) 



and at NNLO: 



■qNNLO p 1 D 1 p 1 p 

-Kll — -Kll,0 + -n-11,1 + -"-11,2 + -Kll,3 , 

r>NNLO _ p i p 

-"■12 ~~ -n-12,1 + -n-12,2 , 

R 22 NLO — -^22,0 + -^22,1 • (2.30) 
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2.3 P, D, F 2 3 and G| partial waves except Df 

For the P and .D waves, the NLO KSW amplitudes [J3] just contain one pion exchange (OPE) and 
at NNLO [T4^ l2~K] they only include in addition the reducible part of the twice iterated OPE. The 
physics behind this is then quite limited and will show up in the phenomenology which, on the 
other hand, is of the same quality as that of the LO Weinberg's scheme results, ref . |T2] . Indeed, at 
LO the potential within the Weinberg's scheme just contains OPE (our order p° contribution) while 
the reducible part of the twice iterated OPE is generated by solving the Schrodinger equation. For 
these partial waves, which are free from enhancements due to unnatural scattering lengths, the 
formalism simplifies since the scaling of the counterterms is just given by dimensional analysis and 
there are no KSW amplitudes of order p -1 . For the elastic ones, 71 = R +Ri+O(p 2 ) and g = O(p ) 
and for the coupled channel partial waves, namely P 2 —F 2 and D\ — G\, Rij = Rij t o + Rij,i + 0(p 2 ) 
and <7i(s) = O(p ), since all the channels start to contribute at the same order p° and as stated 
they are free of unnatural scattering lengths. 

For the partial waves without coupled channels, after performing the appropriate KSW expan- 
sion and matching with the inverse of the KSW amplitude, as done previously, one can write: 



Ra 



Ac 



1 I A Mv 



:P M ^ 



so that, 



,NLO 

u nnlo = Rq + Ri _ (2.32) 



n NLU = Ro, 

i NNLO 



For the two coupled channel partial waves, the chiral expansion is performed in coupled channels. 
Taking into account that R^ = ir^o + R%j,i + 0(p 2 ), as discussed, one has: 



Rn- 



o 



47r(47rAii,o - Mu 2 A\ 2Q + Mu 2 A llfi A 22fi ) 
(4vr + Mz/ 1 A 11> o)(4tt + Mu 2 A 22fi ) - M 2 v x v 2 A\ 



12,0 



#11,1 = t -2 {4ivr(4mA llil (4vr + Mz/ 2 A 22 .o) i 

[M 2 v x v 2 A\ 2f) - (4tt + Mv x A llfi )(4x + Mu 2 A 22fi ) 



R 



+ M(M 2 vlpA\ 2fl - 8iis 2 irA 12A A 12fi (47t + Mu 2 A 22fi ) + pA 2 llfi (4ir + Mu 2 A 2V 
- 2A 2 12fi (-2iMu 2 7iA 22il + p(4n(-2n + Mv 2 A Ufi ) + M 2 u 2 A Ufi A 22fi ))))} , 

16n 2 A 12fi 



12,0 



(4tt + Mz/ 1 A 11i o)(4tt + Mv 2 A 22fl ) - M 2 v x v 2 A\ 



12,0 



#12,1 = {16tt 2 (A 12i1 {M 2 v x v 2 A 

[hPvxVzA^ - (4vr + M^A llfi )(An + Mv 2 A 22 ,o 



2 

12,0 



+ (4vr + Mi^Aii )0 )(47r + Mu 2 A 22fi )) + iMA 12> o(-4pirA llfi + M{v x + v 2 )pA\ 2Q + Mv 2 -nA 22 ± 
+ iMv x v 2 A X \s)A 2%1 - 4pnA 22fi - MvipA n ,oA 22 , - Mv 2 pA 11;Q A 22t0 
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-^22,0 
-^22,1 



+ zz/ 1 A 22i1 (4tt + Mu 2 A 22fi )))} , 

Anj-M^Aj^ + (4tt + Mu 1 Au fi )A 2 2fl) 
(4tt + Mz/ 1 A 11i o)(4tt + Mv 2 A 22fi ) - M 2 v x v 2 A\ 2Q 

1 



2 {4OT(-8iMz/i7r(47r + Mi/iAii )Ai 2 iA 

[M 2 z/ 1 z/ 2 A2 2>0 -(4tt + Mi/iA 1 i |0 )(47T + Mz/ 2 A 22i0 )] 

+ Af 3 i/jfpAi 2)0 + (4tt + Mz/ 1 v4 lli0 ) 2 (4i7Tv4 22il + MpA 2 22 ) - 2MA 2 120 (-2iMv 2 7rA 11A 

+ p(-87T 2 + Mz/ 1 (47r + Mz/ 1 A 11 ,o)A 22i0 )))} , (2 

where the subscript 1 before the period always refers to the channel with lower orbital angular 
momentum and 2 to the highest one. 



Working at NLO implies taking: 



R NLO 



R 
R 



NLO 
12 

NLO 
22 



Rn,o , 

Rl2,0 , 
-^22,0 j 



(2.34) 



and at NNLO: 



tdNNLO 
K u 

tdNNLO 
K 12 

tjNNLO 
rl 22 



Rn,a + Rn,i ) 

Rl2,0 + Rl2,l , 
-^22,0 + -^22,1 • 



(2.35) 



As a check for all the previous expressions for the interacting kernel TZ, we have explicitly verified 
that they are real in the physical nucleon-nucleon region, p real and positive, and furthermore they 
also have the correct KSW order, as they should. 



3 Results and discussion 

In this section we consider the phenomenological applications of the previous scheme to the S and 
higher partial waves. Of particular relevance is to study the triplet S-wave channel, Sf, and its 
mixing with the Df , since here the KSW amplitudes do not converge for p > 100 MeV although 
pions are explicitly included. We first discuss the Sq channel and then consider the Sf coupled 
with the D\ partial wave. After that we turn to discuss the P, D, F 2 and Gg waves and compare 
with other approaches. 



3.1 Sq channel 

We follow the notation of ref.^1] for the KSW amplitudes, where only those combinations of 
counterterms that appear in a given amplitude are shown and are denoted by £j. At NLO we have 
two of such counterterms, £i and £ 2 together with 7 that already appears at LO. We express these 
two counterterms in terms of 7 by performing the ERE in eq. ()2.5j) , with TZ NLO given in eq. (|2.19J) , 
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reproducing the physical values of the scattering length, a s , and effective range, r . We then have: 



-g 2 A M 2 (Q-f 2 - 87771^ + 3ml) (-1 + a s v f + W 2 ™ 2 ^ - 8a s u + a 2 Mis 2 7rr ) 
" ~ 96/ 2 m2(-l + a s v) 2 it 2 ' 

M(^ 7 M (7 - 2m 7 >(-l + a s v) + 8f 2 (-v 2 + 7 2 (-l + o B i/) + 7 z/(-l + q B z/))7r) 
^ 2 ~ 32f 2 mlu(-l + a s u)rt 2 ' l j 

At NNLO there are three more counterterms, £3, £4 and £5. Once the renormalization group equa- 
tions are solved within the KSW perturbative scheme, £5 results to be a higher order counterterm 
and must be set equal to zero at NNLO. We will show below, when discussing the results for 
the Sf — D\ channel after eq. ()3.15|) . that indeed, £ 5 must be left as a free parameter within our 
scheme. Nevertheless, we found after performing fits with £5 free, that it turns out to be negligibly 
small in any case, so that in the following we also make £5 = for this channel as in the pure 
KSW treatment [Hj. We then fix the counterterms £3 and £4 in terms of £1, £2, 7, o, s and tq by 
performing the ERE. The expressions of £3 and £4 in terms of £1, £2, 7, a s and ro are: 

£3 = 77^4 2 I : - T {48/ 2 ^Mz/(-l + a s z/)7r(2 7 3 M + 7 2 M(-4m. + z/) + 8m^vr£2 

1536/ 4 m£z^(— 1 + agUjir 6 L 

- 2 7 wvz/(M + 4m^7r£ 2 )) - 384/ 4 tt 2 ( 7 2 Mz/(1 - a s u) + 7 3 (M - a s Mu) - ~/u(-l + a s v) 
x [Mv - 8m 2 vr£ 2 ) + v 2 [Mv + 4m 2 (-1 + a s v)n&)) + g A M 3 u 2 {-l + a s z/)(6 7 3 - 21 7 2 m^ 

- 7 m 2 (-18 + log(4096)) + m 3 r log(4096))} , 

£4 = ~ qn79f4 rj-T7 ^{384/ 4 m 2 ^(87(-l+a^) 2 (-l + 2vr 2 £ 1 ) + K-4 + 87r 2 £ 1 

+ a 2 z/ 2 vr(-Mr + 8tt£i) + 8a s u(l - 2tt 2 £!))) + 32f 2 g 2 A M(-l + a s i/) 2 vr(12 7 3 M + 7 2 M(-16m^ + 6v) 

+ 2 7 m^(M(3m^ - Av) + 12m Tr z/7r(^ 1 - 2£ 2 )) + m\v(3M + 8m^vr(-3£ 1 + 4£ 2 ))) 

+ ^M 3 z/(-l + a s z/) 2 (48 7 3 - 135 7 2 m^ + 2m 3 (-13 + 41og(16))-4 7 m 2 (-27 + log(4096)))} . (3.2) 

The values for the parameters that we take are / = 130.67 MeV, g& = 1.267, m n = 138 MeV, 
M = 939 MeV. Specifically for the Sq channel, the ERE parameters are a s = —23.714 fm and 
r = 2.73 fm. On the other hand, since a s is so large in this channel, and 7 is around l/a s both 
in KSW as in our approach,* 5 in the following we take directly 7 = 0. 

In the KSW approach at NLO one has the free parameters 7, £1 and £ 2 . In ref. jHj 7 is fixed 
by requiring the presence of a pole in the unphysical sheet in the position required by the ERE, 
just below threshold. Then £ 2 is fixed in terms of £1 to avoid spurious poles that appear at NLO 
(double poles). As a result only £ x is taken as free in the fit to the elastic Sq phase shifts in ref . [T3j . 
The resulting values are: 

Y = -7.88 far 1 , £x = 0.22 , £* = 0.03 , (3.3) 

where those parameters marked with a star are not taken as free ones in the fit. The corresponding 
curve is the dotted one in the left panel of figffl In the approach that we present here we have 



# 5 For example, just take TZ = Ro + R\ in ea. (|2.5() and then perform the ERE, with Rq and R\ given in 
eq.JHH§. 
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one more parameter at NLO, v. Since, as discussed above, we have fixed 7 = and £1 and £ 2 
are calculated from the ERE, eq. fj3.ljl . only v would be free in the fit. Instead of performing 
a fit, and in order to show as well the sensitivity of our results at NLO under changes of the 
subtraction constant u, we show in the left panel of figffl two curves corresponding to v — 200 
(solid line) and v = 900 MeV (dashed line). The solid curve reproduces already very well the 
data for p < 150 MeV. The calculated counterterms are = 0.22(0.25), & = 0.03(0.03) for 
v = 200(900). Taking into account table^we see that the considered variation of v implies a large 
change of the hypothetical cut-off from around 300 MeV to values around 2 GeV. It is also worth 
mentioning the constancy of the values of the counterterms £1 and £2. This is just a consequence 
of the fact that va s » 1 as can be seen from eq. ()3.1j) . 




50 100 150 200 250 300 50 100 150 200 250 300 



p (MeV) p (MeV) 

Figure 1: Phase shifts for the Sq channel. The dotted lines in the left and right panels are the 
NLO/NNLO KSW results respectively. Left panel: The solid and dashed lines represent the NLO 
results of our approach with v = 200 MeV and 900 MeV, respectively. Right panel: The solid line is 
the fit of eq. (|3.5j) . The dashed and short-dashed lines are our NNLO results with v = 500 MeV and 200 
MeV, in order. The data are from the Nijmegen partial wave analysis, ref. |24| . 

Let us consider now the NNLO results both from the pure KSW approach and ours. At this 
order three new counterterms appear, £ 3 , £ 4 and £ 5 . Because £ 5 is a higher order counterterm 
is fixed at zero in ref.^l]. On the other hand, by imposing, as in the NLO case, the absence of 
double and triple poles one can express £3 in terms of £4. As a result only two free parameters are 
left to fit the data. The resulting values from ref. ^1] are: 

7 * = -7.88 fin" 1 , Ci = 0-078 , = 0.03 , = 0.18 , £ 4 = 0.25 . (3.4) 
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The generated curve is the dotted one in the right panel of figQ which indeed reproduces the data 
rather accurately. We now come to our approach. As stated, £3 and £4 are fixed in terms of 7, 
£1, £2, cl s and r by eqs. (J3.2j) . As in the NLO case £1 and £2 are fixed by eq. ()3.1j) and 7 = 0. As 
a result we have only one free parameter, v. After performing the fit to the phase shifts one has 
the values: 

7* = far 1 , £* = 0.25 , £* = 0.03 , £* = 0.21 , £* = 0.24 v = 870 MeV , (3.5) 

corresponding to the solid line in the right panel of figffl which is quite similar to the NNLO KSW 
one with two free parameters although not so close to all the data points. In addition we show 
by the dashed and short-dashed lines those curves that result from our formalism keeping v fixed 
at the values 500 MeV and 200 MeV, respectively. The variation in the results from changes of 
v is similar to that obtained at the NLO, although the quality of the reproduction of data have 
improved. 

We now consider as well the possibility of fixing v and then fitting £1 and £2 while £3 and £4 
are calculated once more from eqs. ()3.2|) . We take the set of values v = 200, 500, 700 and 900 
MeV that are shown in the left panel of figOby the short dashed, dot-dashed, solid and dashed 
lines, respectively, although cannot be distinguished. In addition we also show the dotted line 
corresponding to the NNLO results within the KSW approach from ref . ^1] . We see that basically 
there is no dependence on v although the values of the fitted counterterms change substantially 
and they are in order, for v = 200, 500, 700, 900 MeV: £1 = -0.53, 0.14, 0.56, 0.87 and £ 2 = -0.21, 
0.04, 0.23, 0.36. On the other hand, the reproduction of data is very good. 

Let us now go to higher energies and show results up to laboratory kinetic energies, T i(lfe =350 
MeV (the threshold for pion production is at TJ {, = 280 MeV). This is given in the right panel of 
figCD The solid line corresponds to the solid one of the left panel of the same figure, that is, when 
fixing v = 700 MeV and £1 and £2 taken as free parameters (as shown in the left panel of figOthe 
change from one value of v to another is almost negligible). In addition, the dashed line is the 
solid one of the right panel of figO eq. (|3.5|) . The short dashed line is calculated at NLO within 
our approach with v = 200 MeV (solid line in the left panel of figQ) and the NLO KSW results 
of ref. [Hj are given by the dotted line. 

3.2 Sf — D\ coupled channels 

We treat the counterterms in analogous lines as described above for the Sq channel, although now 
7 and £5 are taken as free parameters. In addition, at NNLO in KSW there is a new parameter 
in this channel, £6, related to the Sf — Df mixing. One should take into account that although 
we use the same names for 7 and the £'s counterterms as in the Sq case, they are indeed different 
[nil E]. As before £1 and £2 are given at NLO in terms of 7, a s and from the ERE. At NNLO 
we then calculate £ 3 and £ 4 from £ 1( £ 2) 7, a s and r , independently of whether £1 and £ 2 are 
either taken as free parameters or fixed at NLO from ERE. The expressions for £1, £2, £3 and £4 
determined from the ERE are the following: 

& = axf2 2( ri ^{-^M 2 (67 2 - 8 7 m^ + 3m2)(-l+a s z/ 1 ) 2 + 12/V(4-8a s z/ 1 

+ a 2 Mz/ 2 7rr )} , 



16 






150 
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(MeV) 



Figure 2: Phase shifts for the Sq channel. The dotted lines are the NNLO results from the pure KSW 
approach ^3]. The rest of the curves are calculated within our scheme. Left panel: The different curves 
are hardly distinguishable. The short-dashed, dot-dashed, solid and dashed curves are our NNLO with 
v taken as 200, 500, 700, and 900 MeV, respectively. £i and £2 are taken as free parameters. The right 
panel corresponds to higher energies with the phase shifts as function of ?} a &. The solid line is the solid 
one of the left panel. The dashed line is the fit given in ea. (|3.5|) at NNLO and the short-dashed one is 
the solid line of the left panel of figH calculated at NLO with v fixed at 200 MeV. The data correspond 
to the Nijmegen partial wave analysis, ref. 24 . 

6 = „ 9 . 2 2 — 2 {M{g\ 1 M{ 1 -2m^)v x {-\ + a s v x ) + %3\-v\ 

32/ 2 m^z/i(-l + a s isi)n 2 

+ 7 2 (-l + a s z/ 1 )+7z/ 1 (-l + a s z/ 1 ))vr)} , (3.6) 



- 2 1 m 1 ,v 1 {M + 4m 7r 7rf 2 )) - 640/ 4 7r 2 ( 7 2 Mz/ 1 (l - a,^) + j 3 (M - a s Mv\) - 7^i(-l + a s v x ) 
x [Mv x - 8m 2 vr6) + v\{Mv\ + 4m 2 (-l + a s v x )ir£ 2 )) + g\M 3 i^{-l + a s ^i)(107 3 - 95 7 2 m w 

- 57m 2 (-15 + log(4096)) + 2ml (-8 + log(67108864))) } , 

6 = - 107520/^M + a, Vl )V {""OAfrPM-l + + 2n%) + , l( -4 + 8^6 

+ a 2 s is 2 ir(-Mr + 8vr^) + 8a 8 z^(l - 2tt 2 ^))) + ll20f 2 g 2 A M(-l + a s z/i) 2 7r(12 7 3 M + 7 2 M(-16m. 
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+ 6i/i) + 2 imn (M(3m w - 4v x ) + Um^n^ - 2£ 2 )) + m>i(3M + 8m 7r 7r(-36 + 4£ 2 ))) 

+ ^M 3 z/x(-l + a s z/i) 2 (16807 3 - 33257 2 m^ - 6tt4(913 + 1248 log(2) + 280 log(4) - 336 log(16) 

- 140 log(256)) - 8407m 2 (-12 + log(256)))} . (3.7) 



At NLO the free parameters are 7, £1, £2, ^1 and v 2 while in the pure perturbative treatment of 
ref.[13[fl4] the latter two are absent. In ref.^3] 7 is fixed by requiring the presence of the deuteron 
pole in the physical sheet and £2 is given in terms of £1 in order to avoid double poles. As a result 
only one free parameter remains at this order in the KSW scheme, £1. This is fitted to the low 
energy Sf elastic phase shifts for p < 80 MeV in ref . J3] with the resulting values ^3] : 

7* = 0.23 far 1 , £ x = 0.327 , £* = -0.0936 , (3.8) 

where as usual the stars indicate that 7 and £2 are not free parameters in the fit. The results are 
given by the dotted line in figOJboth for the Sf elastic phase shifts and for the t\ mixing angle 
which is defined such that S\\ = e s ? cos2ei. The Df elastic phase shifts are presented in the 
figlHland will be discussed below together with the rest of D-waves. The agreement is remarkably 
good for the Sf phase shifts and promising for the e% parameter, in the sense that the NNLO 
contributions are expected to improve the agreement with the e\ data. In the same figure we also 
present the resulting curves from the scheme presented in this work at NLO. So as to reduce the 
number of our free parameters as much as possible we impose that z^ = Z7, that is quite natural 
if we think of the z/j as coming from a cut-off as discussed above. This constraint together with 
eqs. f!3.6|) . which fix £1 and £2 in terms of the ERE parameters, a s = 5.425 fm and r = 1.749 
fm, implies that only 7 and v\ remain as free parameters. We then obtain after the fit to the 
scattering data: 

7 = 0.41 fm -1 , £* = 0.33 , ££ = -0.06 , v\ = v x = 670 MeV . (3.9) 

The resulting curves are the solid ones in figO These curves lie in general somewhat closer to data 
than those of the pure KSW treatment of ref.^3] at NLO, particularly for cm three-momentum 
above 150 MeV. In addition we also consider the sensitivity of our results under a change of the 
subtraction constants z/ 2 and u%. The short-dashed curves, the one lying highest in the left panel, 
correspond to fixing v 2 — V\ = 800 MeV and then 7 is fitted with the result 7 = 0.41 fm -1 with 
£1 = 0.32 and ££ = —0.02. Analogously, the lowest lying curve in the Sf phase shifts, the dashed 
ones, corresponds to taking z/ 2 = v\ = 200 MeV and then 7 = 0.42 fm -1 , ££ = 0.44 and ££ = 0.06. 
Taking into account the lists of values for v given in tabled as a function of an hypothetical cut-off, 
it is clear that a variation of the v's from 200 MeV to 900 MeV, with the best fit at v 2 = v\ = 670 
MeV, can be recast as a large variation of Q, from around 300 MeV up to 1.6 GeV. It is also 
remarkable the constancy of the value of 7 around 0.41 far 1 . 

Now we come to the NNLO order results both at the perturbative level of ref. [13] as well as 
from our scheme. As before, we consider first the pure KSW analysis ^3]- The value of 7 is 
the same as the one at NLO eq. (|3.8|) . since the pole in the perturbative treatment comes only 
from A-i, the leading contribution. Similarly as at NLO, the number of free counterterms can 
be reduced by requiring the absence of double and triple poles that appear in the NNLO KSW 
amplitudes. As a result, £ 2 and £ 3 are expressed in terms of £1 and £4. Thus the free counterterms 
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Figure 3: Phase shifts for the Sf channel and mixing angle e%. The dotted lines are the NLO KSW 
results ^3] . The solid lines are the results from our approach corresponding to the fit of eq. (|3.9|) . In the 
short-dashed lines the subtraction constants v\ and V2 are fixed at 800 MeV. The dashed lines correspond 
to fixing v\ and vi at 200 MeV. The data come from the Nijmegen partial wave analysis, ref . |24j . 

present in ref. [Hj are finally £i, £4, £5 and ^ which are fitted to the S-wave scattering data, at 
this order the Df elastic phase shifts are free of any NNLO counterterm. The counterterm £q does 
not affect the elastic Sf phase shifts and only influences the e\ mixing angle. The dotted lines in 
fig |U correspond to the pure KSW treatment with the set of values [14]: 

7* = 0.23 far 1 , Ci = 0-432 , = -0.0818 , & = 0.165 , £4 = 0.399 , £ 5 = 0.26 , (3.10) 

by fitting the elastic Sf phase shifts. The counterterms with a star are not taken as free parameters 
in the fit as usual. The resulting phase shifts are shown in the left panel of fig0]by the dotted line. 
In the second entry of ref. [T4j only the mixing angle e\ was considered and a fit was performed 
exclusively to the e% data from ref. [21], without considering the Sf phase shifts. The new set of 
values are: 

7* = 0.23 far 1 , & = 0.235 , £* = -0.104 , £ 6 = 0.385 . (3.11) 

The resulting Sf phase shifts are much worse than those from the set of eq. (|3.1()jl . Thus we only 
show the curves with the values of eq. (|3.1()J) and then we determine & by performing a fit to the e\ 
data with the rest of counterterms fixed at the values of eq. ([3.10|) . We obtain then £ 6 = 0.50 and 
the resulting curve for e\ is the dotted line shown in the right panel of figE] As it is well known 
from the results of ref.[Hj, the NNLO results are worse than those at NLO already for p > 100 
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MeV, despite the pion fields being explicitly included in the effective field theory. The results 
shown in ref. ^1] from the set of values of eq. 1)3.11)1 give rise to the same divergent behavior as 
that indicated by the dotted lines of figlU It was also noted in ref . |T3] that this bad behavior was 
due to large corrections from the twice iterated one pion exchange diagrams which are enhanced 
by large numerical factors. In our power counting the input kernel, TZ, is infinitely iterated and 
with it the pion exchange as any other contribution. 

We now consider the results we obtain from our novel non-perturbative approach at NNLO 
order in the expansion of TZ. We follow the same treatment for the counterterms as explained in 
the Sq case, so that £1, £2 and £3, £4 are fixed in terms of the ERE parameters at NLO eq. (jH.fi)) 
and at NNLO eq. (j3.7j) . respectively. Performing a fit to the scattering data, Sf phase shifts and 
ei mixing angle, up to p < 300 MeV, the values we obtain are: 



7 = 0.37fm- 1 , £ = 0.44, £ 2 * = 0.01 , £* = 0.05 , £ 4 * = 0.04 , £ 5 = 0.19 , £ 6 = 0.58 , 

v x = 190 MeV , u 2 = 620 MeV . 

(3.12) 

The results of this fit are presented in the left and right panels of figEJby the solid lines. As we 
see the agreement with data is remarkably good. We now show the sensitivity of the results under 
changes of the subtraction constants V\ and v^- For that, we first impose the constraint z/ 2 = V\ 
as in the NLO case, and perform the fit. The fitted parameters are: 

7 = 0.41 far 1 , £* = 0.34 , £* = 0.00 , £ = 0.06 , £ = 0.16 , £ 5 = 0.25 , £ 6 = 0.51 , 

v x = 500 MeV , v\ = 500 MeV . 

(3.13) 

the results are shown in figE] by the dashed lines. They are quite similar to those of the fit of 
eq. ()3.12)) and essentially identical for the 6\ mixing angle. Let us note that v\ has changed very 
appreciably with respect to eq. ()3.12)) and the fit continues being rather acceptable. Finally, in 
order to have as well a large variation of z/ 2 from the value given in eq. ()3.12)) . we fix v\ = 600 MeV 
and V2 = 200 MeV and fit 7, £5, £6. The resulting values are: 

7 = 0.55 ml- 1 , £ = 0.33, £ = 0.10, £ = -0.19, £ = -0.11, £5 = 0.20 , £ 6 = 0.57 

14 = 600 MeV , v\ = 200 MeV . 

(3-14) 

The results correspond to the short dashed lines. The Sf phase shifts are very similar to those of 
the best fit of eq. ()3.12)) although the e% mixing angle curve is more different than that obtained from 
the fit with z/ 2 = v\, eq. ()3.13)) . Thus, we see that the dependence on the subtraction constants 
v\ and z/ 2 is rather mild and can be reabsorbed to a large extend in new values of the KSW 
counterterms, 7 and £j. 

In the KSW approach, when solving the renormalization group equations one has [T4] . 

mlM 

t* = P-t~ (3 - 15) 

and since p is order p°, given by A^at, then £ 5 is formally a quantity of order p 2 . Hence, within the 
KSW scheme its contributions to the scattering amplitude are order p 2 and start at N 3 LO. This 
situation is the standard one in ref. ^1] and the results from the KSW perturbative treatment with 
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Figure 4: Phase shifts for the Sf channel and mixing angle t\. The dotted lines are the NNLO KSW 
results of ref.jll] with ^5 / and the dashed-dotted line of the left panel when £5 = 0. The solid 
lines correspond to the fit of eq. ()3.12l) . The dashed ones are the fit of eq. Q3.13J) . The short dashed lines 
correspond to the results of eq. (|3,14|) . For further details see the text. The data are from the Nijmegen 
partial wave analysis, ref. |24| . 

£5 = 0, with the rest of counterterms £1 — £ 4 given in eq. (j3.11j) . correspond to the dashed-dotted 
line in figEJ the one that lies above all the other ones in the left panel. We see then that the elastic 
Sf phase shifts diverge for cm three-momenta much smaller than for the £5 7^ case but similar 
to those cm three-momenta where the divergence starts for the e\ mixing angle, which indeed is 
independent of £5 at NNLO in the KSW approach. 

Now, let us discuss why in our approach we should take £5 as a free parameter without taking 
into account eq. (j3.15j) . The point is the following. The order of the counterterms in KSW, in 
particular those given rise to £5, are determined by comparing with the effective range expansion in 
the pionless effective field theory once the power divergence subtraction (PDS) scheme is adopted 
[T3] . This comparison is exact, even if the series are truncated, since both of them implies an 
expansion in powers of the cm three-momentum and the comparison is performed order by order. 
In this process one formally books the scattering length as order p~ x and the shape parameters 
r n as order zero, proportional to 1/Ajvat- Once this is performed, those operators that enter in 
the Lagrangian up to some order can be determined, and then the corresponding KSW partial 
wave amplitudes up to the same order can be calculated as well. As we have already explained, 
these amplitudes are taken as the input of our approach in order to fix TZ by the matching process 
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Figure 5: Phase shifts for the Sf channel and mixing angle t\. The dotted lines are the NNLO KSW 
result ^1] with £5 7^ 0. The short dashed lines represent the NLO results from eq. (|3.9j) . The dashed lines 
are the NNLO fit to data for p < 300 MeV given in ea. (|3.12|) . The solid lines are the new fit of ea. (|3.16j) 
for p < 400 MeV. The data correspond to the Nijmegen partial wave analysis, ref . |24j . 

outlined above in sec|21 In the pure KSW treatment one requires the amplitudes being independent 
under changes of the subtraction point so that the renormalization group equations (RGE) follow 
for the different local counterterms. In practice, even when pions are removed and only local 
operators remain, this guarantees that the scattering partial wave is subtraction point independent 
only perturbatively, in the sense that this only happens when the Schrodinger equation is solved 
perturbatively.* 6 But this is precisely the point we want to avoid in our formalism so as to take 
care properly of the large 2M/p factors that appear from the two nucleon intermediate states in the 
unitarity bubbles. At this point our formalism is an hybrid between that of Weinberg j^j and the 
KSW formalism [13J. On the one hand, we perform an expansion of an interacting kernel 1Z as in 
ref. jS] but on the other we do that by considering directly the scattering amplitudes as in ref. [TB] . 
Indeed, there is a residual dependence in our amplitudes on the subtraction constant u, expected 
to correspond to higher order operators in the KSW Lagrangian, similarly as in the Weinberg 
approach where there is cut-off dependence that is expected to become softer as the order of the 
calculation increases, as also occurs in our case as discussed above. To sum up, our scheme makes 
use of the KSW power counting in order to obtain the Lagrangian, employs Feynman diagrams 

# 6 This acquires a clear meaning if one thinks of the pionless effective field theory, where to solve the Schrodinger 
equation is a trivial task. 
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to calculate from this Lagrangian the KSW amplitudes at some order (where the order of any 
combination of counterterms is determined directly from the ones of the original counterterms in 
the Lagrangian, without making use of perturbative RGE arguments). Then, these amplitudes 
are matched with the general expression of eq. (|2.5|) so that the interacting kernel TZ is determined. 
Thus, as a result, one ends with partial wave amplitudes with the unitarity bubbles resummed as 
in Weinberg's scheme but in an analytical way based on unitarity and analyticity. 

Encouraged by the rather good fits obtained from our approach up to p < 300 MeV as depicted 
in figsIBl and 01 we now show in figEl our results for higher energies, namely for T} a b < 350 MeV, 
as similarly shown in the studies of nucleon-nucleon scattering within the Weinberg approach in 
ref. |T21 . Let us note that for p = 300 MeV one has Ti ab ~ 190 MeV. We present in fig[5] several 
lines. The dashed lines correspond to the solid ones of figlH with the parameters given in eq. (j3.12|) . 
We see, as already shown in figJH that the agreement with data for T[ ab < 190 MeV is very good 
but from TJ a & = 200 MeV starts deviating from data. The dotted lines are the NNLO KSW results 
from ref. [14 with £ 5 ^ 0, already presented in fig|4]with the same type of line. We also show with 
the short-dashed lines our NLO results given in eq. (J3.9|) with only two free parameters in the fit. 
We see that they follow closely the trend of the Sf phase shifts in all the energy interval although 
the agreement is not so good for the e\, similarly as also happens in the KSW treatment at NLO. 
Finally by the solid line we present a new fit up to p < 400 MeV, where £i and £2 are fitted, 
instead of begin fixed by the ERE at NLO, eqs. (j3.fi)) . On the other hand, £3 and £4 continue to 
be given in terms of £1, £2 and 7 from eqs. ()3.7|) . The resulting parameters are: 

7 = 0.51fm- 1 , 6 = 0.27, 6 = 0.06, £ = -0.10, £4 = -0.05 , £ 5 = 0.17 , £ 6 = 0.46 , 

v x = 180 MeV , v 2 = 560 MeV . 

(3.16) 

The reproduction of the data is quite good, specially for the mixing angle e\. Nevertheless, we 
observe some deviation from data in the low energy region around 7} a & ~ 100 MeV (p ~ 215 
MeV) which does not appear in the other fits obtained by fitting data only up to p = 300 MeV or 
Tiab ^$ 190 MeV, that is for lower energies. This can be interpreted as an indication that we are 
forcing too much the capability of our approach accordingly with an estimated Ajvtv around 400 
MeV. This is also supported by the divergence of the fit of eq. ()3.12|) for Ti ab above 200 MeV, solid 
lines in figEJand dashed ones in figH3 as well as from the results shown in fig|Ulwhich are discussed 
in the next paragraph. In this sense the fit given in eq. (j3.16|) is just a way to fit data up to rather 
high energies. Nonetheless, it is reassuring that the values of the parameters, although there are 
now two free parameters more, are very similar to those obtained above for lower energies, e.g. in 

eq.flTH). 

We also show in figElthe absolute values \1Z NL °\ and \jz NNLO \ for the Sq channel (solid lines) 
and the same, although divided by four to keep the lines on the same scale, for the matrix element 
I -Rnl of the Sf channel. The thick lines refer to the NNLO calculations and the thin ones to NLO. 
The convergence properties are quite good in a broad range of the center of mass three momentum 
p, with strong divergences for p around A^tv — 400 MeV. This scale is the one expected for the 
KSW EFT [13] . although in the end this EFT does not converge in the triplet channels for p > 100 
MeV, as already discussed. We see that within our scheme, at the same time that we keep the 
KSW power counting, we are able to fulfill these expectations. This is the main aim of the present 
investigation. 
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Figure 6: \TZ NLO \ (thin solid line), \K NNLO \ (thick solid line) for the S£ channel and \B$f° /4| (thin 
dashed line) and \R\i NLO /4\ (thick dashed line) for the Sf channel, as a function of the cm three- 
momentum. 



3.3 P, D, Ff and G% partial waves 

Up to NNLO in the KSW approach, order p, the physics included in the description of partial 
waves higher than S waves, except for the mixing between the Sf — D\ ones, consists only of the 
one pion exchange, order p°, and the reducible part of the twice iterated one pion exchange, order 
p. While the former contribution constitutes the LO one of the Weinberg's scheme (HI the 
latter is generated by solving the Lippmann-Schwinger equation. Hence, the NNLO results in the 
KSW approach are included in the leading ones of Weinberg's scheme. It is clear then that we are 
neglecting important contributions for the higher partial waves, particularly for the P and D waves, 
since from the results of the Weinberg's approach we know that two pion exchange irreducible 
diagrams, counterterms for the P- waves at order p 2 , as well as the irN counterterms (which are 
saturated to a large extend by the A isobar) are important contributions that would rise in higher 
order calculation within the KSW, beyond NNLO. Hence a full N 3 LO calculation from the KSW 
Lagrangian should be pursued and then used within our approach to offer a more complete study 
of these higher partial waves. For the F and G waves and mixing parameters €2 and €3, pion 
exchange dominates and the aforementioned extra contributions are not so important, see e.g. |25j. 
This is also clear from our results given in figs|7]and|Hl At the order we are working in these partial 
waves there are no counterterms and the only free parameters are the subtraction constants z/s, 
one for each partial wave. In most of the channels, they take arbitrarily large or negative values, 
that is, any value with modulus typically above A xp t — 0.7 — 1 GeV gives essentially the same 
results. Our curves are indeed quite similar to those obtained in the Weinberg's approach at LO, 
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see [12] ■ Our results at NNLO, solid lines, improve those of the NLO, dashed lines, except for 
the D\ where, though the NNLO is better at low T/ &, they depart from data more than the NLO 
ones for higher energies. They also improve the results from the pure KSW amplitudes, dotted 
lines, although in these expected, the resummation effects are not so important as in the 

S-waves. 

Since most of the subtraction constants v for the analyzed P, D, F and G waves take arbitrarily 
large absolute values we would like to show that in the limit when | — > oo our approach reduces 
to the Inverse Amplitude Method [22] ■ This method would consist of expanding the inverse of a 
partial wave, following in our present problem the KSW power counting, and then calculate the 
partial wave by inverting exactly the previous expansion. For instance, we have given in eq. (j2.15|) 
the expansion of the inverse of the elastic S$ partial wave within the KSW power counting (in 
ref. [22] one uses for meson-meson scattering the standard chiral perturbation theory counting 
and for pion-nucleon processes the chiral perturbation theory counting is supplied with the heavy 
baryon one). Then the inverse amplitude method would imply: 



T, 



A \ , ( A\ - A X A_^ 



A_J \A 2 _J V A ~i 



+ 



-i 



(3.17) 



Its generalization to coupled channels is straightforward by making use of a matrix notation. For 
example in the Sf — Df coupled channels one should just invert exactly the matrix (z^ 5 ^) -1 
whose matrix elements are given in eqs. (j2.25J) . and so on for any other partial waves. This kind 
of results are also usually referred as Pade resummations. 

Let us see that our formalism reduces to the inverse amplitude method when |z/| — > oo. Consider 
first the elastic case. Then, the contribution of order p l in the KSW power counting, from the 
expansion of 1 /1Z + g in eq. ()2.5p . has one piece involving the Ri component of 1Z. Isolating Ri one 
has: 

R l = -R 2 [(l/A) t + g l ] +..., (3.18) 

where gi is now the contribution of the g function of order p l and (1/A)i is that of the inverse 
of the KSW partial wave. Since R scales as 1/z/, see eq. (l2.18j) for the case of large scattering 
lengths and ()2.31|) for the case of natural ones, then Ri has a contribution scaling as 1/z/ 2 . This 
is the only contribution for i = 1 and then R\ scales as 1/z/ 2 . Let us demonstrate by induction 
that each R^ with i > 1 scales as 1/z/ 2 , while Rq scales only as 1/z/. To show that, let us consider 
the dots which indicate other terms from the expansion of 1/TZ involving R m with 1 < m < % — 1. 
We can write these contributions as: 

i— l R ai R a2 R am 

E(-ia\a,2 ...am. ^1 fcgJJJ k m /q i q\ 

Ly k 1 k 2 ...k m ^ai+a 2 + ...+a m -l ' {O.LJ) 
m=l 

where k p (I < p < i — 1) indicates the order of the corresponding factor Rk p , in addition all the 
k p are different from each other and fulfill aiki + 02^2 + ••• + a mk m = %■ Because each of the R^ v 
in eq. ()3.19j) scales at least as 1/z/ 2 one can count easily the dominant power of v from eq. (j3.19j) 
when |z/| — > 00. One simply has: 

. ai+a 2 +...+a m -l i 

- (3 20) 

l/ 2(a 1 +a,2+...+a m ) ly a 1 +a 2 +...+a n +l ' \ ■ J 
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But the sum a\ + a-i + ... + a n + 1 > 3. To show this, let us consider first those terms in the 
sum of eq. ()3.19j) with all the factors the same and equal to R kl . It follows then that a\ > 2 since 
1 < k\ < i — 1. The rest of terms will have at least two different factors, then there must be two 
a's, let us say a\ and 02, different from zero and every one > 1. Then the scaling of Ri in v when 



\v 



— > 00 is dominated the one from eq. ([3.18|) and goes like 1/z/ 2 . 
Let us consider the expansion of 1/TZ + g in eq. (|2.5|) in powers of \ jv in the limit \v\ — > 00 with 
7Z determined up to order p n in the KSW counting. We also write for g only the sum Ylk=o 9k 
since higher order terms are just relativistic corrections and numerically are negligible: 




\ m=l u 



- -I 



2T Lf+ljLf+1 = - [ - + ^g k ) = - [ - - ^ ^ + 0(„-) + ^g k 

k=0 

(3.21) 

Where we neglect the 0{y x ) terms compared to those explicitly shown in the equation above 
which are O(v ) and are the only ones that survive in the limit \u\ — > 00. Thus, as a result of 
the matching process detailed in sec|2]to obtain the Ri functions that we apply now for \u\ — > 00, 
one has between the brackets of eq. (|3.21|) the expansion of the inverse of the KSW amplitude up 
to order p n . So that the resulting T L 2s+i L 2s+i partial wave from eq. ([2.5)1 for \v\ — > 00 is the one 
given by the inverse amplitude method eq. ([3.17|) . For finite z/s this is no longer the case. The 
neglected terms of order v~ x and higher in eq. (j3.21|) guarantee the matching with the expansion 
of the inverse the KSW amplitudes up to order p n and also give rise to higher order terms beyond 
those considered in the matching process. It is only necessary a little thought to extend this 
discussion for the elastic case to the coupled channel one by making use of a matrix language. 



4 Conclusions 

We have established a new analytic expansion in order to treat systematically nucleon-nucleon 
interactions. Analyticity and unitarity properties of the inverse of a partial wave amplitude 
are used to resum the unitarity bubbles in terms of an on-shell interacting kernel 71 and the 
unitarity loop function g. The former is determined by performing an expansion with the KSW 
power counting of the aforementioned general expression for a partial wave and by matching this 
expansion with that of the inverses of the pure KSW amplitudes up to a definite order. In practice 
we have performed this matching with NLO and NNLO KSW amplitudes [T3*l I14j. As a result, 
a hybrid scheme emerges that treats pions non-perturbatively but in harmony with the KSW 
power counting, and that restores the expected range of convergence of the KSW EFT about 
Ajvtv — 400 MeV. It reminds also the Weinberg's approach [5 j in that a perturbative expansion is 
performed for an interacting kernel instead of making it directly to the partial waves. In our case 
this interacting kernel is TZ while in Weinberg's approach is the potential, V. However, it is worth 
stressing that while 7Z is on-shell, V is off-shell and this is an important point that allows our 
formalism to be purely analytic and fairly more simple than that of ref. [5]. It is also important 
to remark that for the S-waves we have achieved a very good agreement with data, both for the 
phase shifts as well as for the mixing angle e\, up to the opening of the NNn threshold. 

Further applications of the present scheme should be pursued, particularly to three body prob- 
lem where the disposal of analytical methods in the two-nucleon sector is well worth [2E| and 
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furthermore a KSW N 3 LO calculation should also be performed within our scheme particularly 
for a more complete study of the P and D partial waves. Finally, another issue is to extend this 
knowledge of the nucleon-nucleon interactions in vacuum to nuclear matter [2H| 130 j. 
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Figure 7: Phase shifts for the P±, Pq, Pf, P| > an d £2 from left to right and top to bottom, respectively. 
The dotted line, when present, is the NNLO KSW result |14j . The dashed line represents the NLO results 
from eq. IJ2.5jl . The solid lines are the NNLO results. The data correspond to the Nijmegen partial wave 
analysis, ref . |24j . 



29 



Ti ab (MeV) 



Tia b (MeV) 




150 200 250 300 




T lab (MeV) 



T lab (MeV) 




T lab (MeV) 




T lab (MeV) 



Figure 8: Phase shifts for the D\, D\, D%, D%, G| and from left to right and top to bottom, respectively. 
The dotted line, when present, is the NNLO KSW result |14j . The dashed line represents the NLO results 
from ea. (|2.5|> . The solid lines are the NNLO results. The data correspond to the Nijmegen partial wave 
analysis, ref . |24j . 
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